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Abstract 

Incorporating the QCD axion and simultaneously satisfying current constraints on the dark 
matter density and isocurvature fluctuations requires non-minimal fine-tuning of inflationary pa- 
rameters or the axion misalignment angle (or both) for Peccei-Quinn symmetry-breaking scales 
fa > 10^2 QgY^ To gauge the de gree of tuning in models with many axion-like fields at similar 
symmetry-breaking scales and masses, as may occur in string theoretic models that include a QCD 
axion, we introduce a figure of merit that measures the fractional volume of allowed parameter 
space: the product of the slow roll parameter e and each of the axion misalignment angles, ^o- Foi' 
a single axion, JF < 10^^^ is needed to avoid conflict with observations. We show that the fine tun- 
ing of becomes exponentially more extreme in the case of numerous axion-like fields. Anthropic 
arguments are insufficient to explain the fine tuning because the bulk of the anthropically allowed 
parameter space is observationally ruled out by limits on the cosmic microwave background isocur- 
vature modes. Therefore, this tuning presents a challenge to the compatibility of string-theoretic 
models with light axions and inflationary cosmology. 
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Introduction 

The axion is today considered a leading solution to the strong-CP problem as well as a 
dark matter candidate and a possible feature of string theory. The axion solution to the 
strong CP problem entails the introduction of a global U{1) symmetry, known as 

the U{l)pQ or Peccei-Quinn (PQ) symmetry that is spontaneously broken to a discrete 



subgroup Z]\f, generating a pseudo-Nambu-Goldstone boson, known as the axion 
Through non-perturbative interactions, the axion, </>, obtains a periodic potential with N 
distinct minima separated by A0 = 27r fa where fa is the PQ symmetry breaking scale. 
In the original proposal, fa was set to the electroweak scale, fa ~ TeV, resulting in an 
axion with large couplings to matter fields; this was quickly ruled out. It was soon realized, 
however, that the symmetry-breaking scale could be much higher, reducing the couplings to 
matter and rendering the axion difficult to detect. This possibility, known as the "invisible 
axion," is the realization addressed in modern discussions js,!?]. 

In the big bang inflationary model, two scenarios of axion production are possible. If the 
symmetry-breaking scale fa is less than the Hubble parameter Hj at the end of inflation, the 
axion is non-uniform when it is generated and the universe breaks up into tiny regions sep- 
arated by domain walls (assuming > 1); the domain wall energy dominates the universe, 
resulting in an unacceptable cosmology. We refer to such axions as low-fa axions. Alterna- 
tively, if PQ symmetry breaking occurs before inflation [high- fa axions), the axion field is 
made uniform throughout the observable universe during infiation and set at some random 
value, typically displaced from the minimum of its potential. The fractional shift from the 
nearest minimum, A0/(27r/a), known as the misalignment angle, determines the amplitude 
of the field once it starts to oscillate about the minimum. The oscillations of the coherent 
Bose- Einstein condensate are pressureless, like dust; this non-thermal energy is a candidate 
for explaining cold dark matter l8|. Numerous experiments have been attempted to detect 

nrin 

the axions produced by this "misalignment mechanism" (see, e.g., [9|, [lO, iJJj ) . Assuming a 



value of the misalignment angle 9q ~ which is typical, a value of < 3 x 10^^ GeV is 

required to avoid the overproduction of axionic dark matter via misalignment Q, [3], and 
a tighter bound of < 3 x 10^*^ GeV arises when the radiation from axionic strings is taken 
into account A lower bound also arises from axion couplings to photons and other 

particles: /„, > 10^ GeV is required to evade other astrophysical and laboratory constraints 
(see, e.g., (uj). 
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It is possible to accommodate a high-/a QCD axion with fa > 10^^ GeV only if we live in 
a rare region of space with Oq exponentially small compared to unity. To explain the unlikely 
value, anthropic selection must be invoked based on the argument that a high axion density 
much greater than the observed dark matter density) is hostile to the development of life 
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In addition to the misalignment, quantum fluctuations induced in the axion field during 
inflation influence axion field oscillations; even if 6'o = 0, inflationary fluctuations displace 
the field from its minimum and lead to the production of some axion particles. The rms 
fluctuation in the axion misalignment angle due to inflationary perturbations is 

Hi 



(1) 



Therefore, to suppress the production of axion particles after inflation, for > 10^^ GeV, it 
is not only necessary to have an unlikely misalignment angle, but also the inflationary energy 
scale must be small. As shown below, this requires finer tuning of the inflationary slow-roll 
parameter, e, than is needed to solve the cosmological flatness and horizon problems. 

The axion scenario may be realized in string theory. As pointed out by many authors 



(e.g., [18|, ll9|, |20|) , compactifications in string theories always generate PQ-like symmetries 
and generically produce a multitude of axions and axion-like fields (ALFs). The number 
and properties of the ALFs depend on the specific model being considered. For the case 
of the heterotic string, there are two classes of ALFs. The so-called "model-independent" 
axion arises from the antisymmetric tensor field -B^jy, where fi, v are the usual space-time 



indices, that plays a key role in anomaly cancellation 



2l|. This axion, which has analogs 



in other asymptotic limits of string theory and does not depend on compactification, has a 
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231]. In the other class of ALFs 



symmetry breaking scale fa of order 10^^ GeV 
are the "model-dependent" axions whose properties depend upon the compactification of 
the manifold. These arise from B^j^ where i, j are internal indices, and typically have /alf 
values between 10^^ and 10^*^ GeV. It is possible to construct models with f alf somewhat 
lower (for example, in warped heterotic string theory; see 2J]), though these models are not 
generic. Cosmological constraints apply equally to the QCD axion and to ALFs with high 
symmetry- breaking scales /alf ^ 10^^ GeV. 

ALF symmetry-breaking scales /alf and masses (m^LF ~ ^\lfI Ialf, where A^lf is 
a scale set by the ALF coupling to instantons) may span a large range, with each ALF 



3 



associated with its own gauge group [19|. This range can conceivably include ALFs with 
^ALF values ranging from the QCD scale to the string scale. If all ALFs are very massive 
(AyiLF ^ Hj), they begin oscillating before inflation, producing particles that are inflated 
away and are, therefore, harmless. However, if string theory can produce a QCD axion with 
a large symmetry-breaking scale and a very small mass (< 1 eV), it is reasonable to 
expect the spectrum to include other ALFs whose mass and Jalf are within a few orders of 



magnitude of those of the QCD axion. Further, as argued for example in Ref. [18|, [191, |20 |. 
string theory models could produce many ALFs with Jalf values near the GUT scale (> 10^^ 
GeV). The masses of these ALFs should be homogeneously distributed on a log scale, with 
possibly several ALFs per decade of energy [l^. Like the QCD axion, all of these ALFs 
that begin oscillating after the end of inflation (i.e., those with Jalf ^ Hj > A.alf) can 
contribute significant densities of dark matter, and can produce isocurvature perturbations 
observable in the cosmic microwave background (CMB). 

Previous works have made use of the possibly large number of string theoretic ALFs to 
construct models of inflation 25| and quintessence 26|. In those constructions, the effect 
of the ALFs is to contribute to accelerated expansion when the axion field is overdamped 
by Hubble expansion. In this work, we consider the limit in which string theoretic ALFs 
have underdamped oscillations which act as dark matter. We show how the existence in 
string theory of many axion-like fields can exponentially exacerbate the cosmological and 
fine-tuning problems for axions in the context of inflationary cosmology. 

To quantify the fine tuning required to evade cosmological constraints in the presence of 
multiple ALFs, we introduce a figure of merit JF which measures the fractional volume of 
allowed parameter space including the tuning of ALF misalignment angles and the slow-roll 
parameter e of the inflationary model. Both Oq and e are 0{1) in an untuned model. We 
define: 

^ALF 

:F^el[9i (2) 

i=i 

where ualf is the number of ALFs, labeled by j, each having its own misalignment angle 
6q. The slow-roll parameter is 



where mpi ^ 2.4 x 10^^ GeV is the reduced Planck mass. This parameter measures the 
flatness of the inflationary potential V and is related to the inflationary equation of state 
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w = pI p, with e = |(1 + w). 

Obtaining the 60 e-folds of inflation required to solve the flatness and horizon problems 
requires e < 10~^, a modest fine-tuning. The amplitude of density fluctuations produced 
during inflation is 5p/ p ~ Hj / {m-p\^). Since observations require 5p/ p ^ 10"^, we must 
have e = 10^°(-ff//mpi)^. This requires Hi < 10~^mpi; making Hj smaller requires more 
fine tuning of e than is necessary for cosmology. For each order of magnitude reduction in 
Hi, two orders of magnitude additional tuning of e is required. (Some have argued that 
very small Hi is essential to incorporate inflation in string theory; if true, this could be 
interpreted as a sign that large Hi is ruled out or that inflation and string theory are a 
poor fit. However, either interpretation is premature. Small Hi models have proven to be 
problematic and string theoretic examples with large Hi have been constructed 28|. 
So, without a stronger microphysical argument, treating e as a free parameter with uniform 
prior is the conservative choice.) Defining the figure of merit as in equation ([2]) provides 
a measure of the fractional volume of parameter space where each do is independent and 
stochastic and e is also independent of each of the ALFs. 



A sing" 



paper 



e QCD axion corresponds to fine tuning at the level of JF < 10 . In a companion 



29|, we show that, in the majority of the parameter space, the tuning is necessary 
primarily to avoid isocurvature modes. Since there is no known reason why an observable 
isocurvature contribution makes the universe uninhabitable, anthropic reasoning cannot 
explain the required fine tuning. In order for the dark matter density constraint to dominate 
the limit on the misalignment angle, one must carefully tune e, which is also unsupported 
by anthropic considerations. 

In this work we will show that each additional ALF with comparable mass exponentially 
exacerbates the fine tuning JF. The result is then used in a Bayesian analysis to compare 
models in which many ALFs are produced with those in which they are not. We also 
find that, as in the case of the single axion 29||, anthropic reasoning is ineffective because 
most of the anthropically allowed parameter space is ruled out observationally due to the 
isocurvature constraint. 

Cosmological Constraints on Axions and Axion-Like Fields 



The natural value for fa (or Jalf) in string theory is of order the string scale, > 10 



GeV 



15 



191 ]. corresponding to a high-/a axion. The axion field is made uniform and 



frozen at some misalignment angle 6*0 during inflation, and the field remains frozen until 
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the Hubble scale is of order the axion mass, m^, at temperature Tqsc- The precise form 
of the temperature-dependent ma{T) relies on the relationship between Tosc and the scale 
A at which instanton effects become important; two limiting cases for the QCD axion are 



discussed in 



12j. For fa < 0.26 (A/200 MeV)2mpi, T^sc is greater than Aqcd- In this 



regime, the axion mass is temperature dependent because the QCD instanton density is 
temperature dependent - this density determines the mass. For larger /„, Tosc < Aqcd so 
rrialT) is temperature- independent. The temperature dependence of a QCD axion is not 
straightforwardly generalizable to ALFs because it depends upon the details of the gauge 
couplings. When calculating the cosmological abundance of ALFs, we treat the mass as 
being independent of temperature; our conclusions are not sensitive to this assumption. 

The total density of the QCD axion plus any additional ALFs is currently observationally 
constrained by the dark matter density and the fraction of isocurvature perturbations in the 
CMB. These constraints can be used to set limits on the axion's fundamental scale {fa or 
fALp) and its post- inflation misalignment scale Oq, as well as on Hj. Detailed derivations 
of the density and isocurvature constraints can be found in, e.g., ll^- Here we briefly 
describe the origins of the constraints and quote the relevant formulae. Except for the 
differences due to the treatment of the temperature-dependent mass, the formulae apply 
equally well for the QCD axion and for ALFs. 

Following |l4|, we parametrize the axion density by the late-time axion energy density 
per photon, ^a, defined by: 

^ ^ PajTo) _ ma{To) pajTosc) s{To) 

n^(To) niaiTosc) n^{To) s{Tosc) 
where s is the entropy density, pa is the mass density of axions, is the number density of 
photons, and Tq refers to the present day. C,a is related to Qa, the ratio of the axion density 
to the critical density, by Qah"^ ~ 'Ca/26 eV with h the dimensionless Hubble parameter 
{h = -ffo/100 km/s/Mpc). The observational constraint on dark matter ^cdm ~ 2.9 eV 
places an upper limit on the axion density. 

Inflationary fluctuations in the axion field that contribute to the axion particle density also 
lead to isocurvature perturbations in the CMB. Limits on the fraction of CMB perturbations 
that are isocurvature can, therefore, place constraints on an axion field existing during 
inflation. To test against this limit, we calculate the ratio of the average power in the 
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isocurvature component to the average total power in CMB temperature fluctuations: 

m/T)u 



an 



(5) 



30|. In terms of 



Current observational constraints from CMB observations limit < 0.072 
axion parameters, the density and isocurvature can be written 

^ A{9l + al)G, (6) 

(Tg(26'o +CTg)G . (7) 



" 25 m/T)l,)) 

Here, A is the instanton scale of the axion field (or ALF), is the average misalignment 
angle within the horizon, and is the matter energy density per photon. We neglect 
dimensionless factors of order unity that correct for anharmonic effects in the axion potential 
and uncertainties in the temperature dependence of the axion mass. The dependence of 
and Ua on (or Jalf) is absorbed into the function G, which accounts for different cases 
of the temperature dependence of the axion mass. We have assumed that the ALF mass is 
temperature-independent near T = Tosc, so Galf ~ 4.4(/alf/'"^pi)^^^- In the case of the 



QCD axion, Gqcd ~ 2.8(Aqcd/200 MeVy/^ifJmp,V/^ 



and 4.4(/a/mpi)3/2 for > 0.26(Aqcd/200 MeV)2mpi 1^, where Aqcd ~ 78 MeV. New 
calculations of the axion mass temperature dependence have recently been published that 



or < 0.26(Aqcd/200 MeV)^mpi 



slightly alter these scalings [13|, but the small quantitative changes do not affect our result. 

In Figure[T]^a), we plot the cosmological constraints on the QCD axion from measurements 
of dark matter density and CMB isocurvature for three choices of the misalignment angle: 
6*0 = 10~^°, 10~^ and 1. As can be seen in equation ([6]), the form of the constraints depends 
strongly upon the relationship between the (unperturbed) post-inflation misalignment angle 
6o and its rms fluctuation amplitude ae due to inflationary perturbations. For the density 
constraint, the presence of the factor (^q + ag) implies that, when ag << 6q, the constraint 
depends only weakly upon Hj and so primarily sets an upper limit on fa- In the plot, this 
produces a horizontal line above which models are ruled out for 6q = 10~^ and 1. However, 
when ag » 6q (as for 6q = 10~^° in the figure), the constraint mainly rules out high values 
of Hj. For the isocurvature constraint, the dependence on the inflationary perturbation 
amplitude is stronger, since 9o and ae appear in the combination o"g(26'g + crj). Therefore, 
the constraint strongly restricts the upper bound on Hj, as can be seen in the nearly vertical 
part of the 9q = 10~^ curve. When 9q is very small (6*0 < 10~^), the isocurvature constraint 



is everywhere more stringent than the density constraint. For these values, the density 
constraint's cut-off of high fa values lies above the plot (as in the Oq = 10^^° contour) and 
its limits on Hj are less stringent than those of the isocurvature constraint. As 6q increases 
toward 0{1), the isocurvature constraint is still the tighter constraint unless the Hubble 
scale of inflation is very low. 
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FIG. 1: (a) Cosmological constraints for tlie standard QCD axion alone for three values of the mis- 
alignment angle ^o- Each solid curve represents the combined constraint from dark matter density 
and CMB isocurvature fraction for 9q = 1, 10"'^, 10^^'^. For each solid curve, the shaded region 
to the right (and/or above) the curve indicates the region ruled out by cosmological constraints. 
The dashed and dotted curves indicate for = 10~^ the constraints from dark matter density 
and CMB isocurvature fraction, respectively, (b) Combined cosmological constraints for the QCD 
axion plus a single ALF with Aalf / Aqcd, and with 9q^^ = 9q^^ = 10~^. For each solid 
curve, the shaded region to the right and above is the region ruled out by a combination of the 
dark matter density and CMB isocurvature constraints. In both (a) and (b), the white region in 
the lower left indicates the parameter space that is allowed for ^ 1- The hatched area in the 
lower right is the low-fa region (fa ^ Hj) not considered in this paper. 

Dark matter density and isocurvature constraints only apply if the axion or ALF does 
not decay away in the early universe. The QCD axion is a long-lived particle. Its decay into 
photons can occur through an interaction term (e.g., see [3l|) in the Lagrangian: Cint = 
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ga.yyaE-'B + ... where ga-yy ~ ajinfa/N) [up to constants], a is the axion field, and E and 
B are the electric and magnetic field vectors, respectively. The decay rate to two photons 
can then be written r^^^ ~ (yf^^^m^/(647r), with a decay time r = h/Ta-y-y- Comparing this to 
the age of the universe, Tq ~ 4.3 x 10^'' s, we find that the condition for an axion to survive 
to the present day is: 



This condition is easily satisfied for QCD axions with the symmetry breaking scales we are 
considering here; their decay into photons can, therefore, be ignored in our calculations. 

For ALFs with high instanton scales A^^^, however, decay into gauge bosons may become 
important. Using equation ([8]) as an estimate leads to an effective upper limit on the 
scale Aalf for an ALF to exist at non-negligible densities today. For the lowest fa values 
considered here (/„ ~ 10^ GeV), the upper limit is Aalf ^ 100 GeV; at fa ~ 10^^ GeV, as 
expected for general string theory ALFs, the limit is Aalf ^ 10*^ GeV. Hence, we assume 
Aalf < 10® GeV in our analysis. 
Additional Axion-like Fields 

We now consider the effect of multiple ALFs with masses comparable to that of the QCD 
axion and fa ^ 10^^ GeV on cosmological constraints. For only a single QCD axion, the free 
parameters are the misalignment angle 6q and the energy scale of inflation. Each additional 
axion has its own 6^ and Aalf- 

Since the constraints are strongly dependent on the misalignment angle of the axion fleld, 
if any one of the ALFs has a large Oq, it dominates the density and isocurvature contributions. 
Similarly, ALFs with larger values of Aalf at the same 6q also dominate. For the purpose 
of illustration, we show in Figure [H^b) the combined constraint (density and isocurvature) 
for a QCD axion plus one additional axion at a range of Aalf values, setting the 9q values 



As Aalf/Aqcd increases, the ALFs are constrained to lower values of fALF and Hf, 
recall that e scales as H]. In the regime in which ae << Oq (i.e., negligible inflationary 
perturbations), the bound on scales as A^^^^q "^^"^ when ^a is held constant or as Aj^/^pO^ 
when aa is held constant. In the opposite case where ae » Oq, the value of Hj must be 





to be the same 



lowered to evade the constraints. The upper limit goes as Hj 
constant and Hj ~ ^vilf-^alf ^^r held constant. 



A 



^fa^^ for ^a held 
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The contributions of axions to the dark matter density and to the isocurvature pertur- 
bations are both additive, so each ALF contributes independently to a tightening of those 
constraints. If the ALFs are produced with a range of A^^^ values, the constraint is domi- 
nated by the subset of ALFs within about an order of magnitude of the largest A^^^. Among 
these, the ALFs with the highest 6q will provide the dominant contribution to the constraints 
at the same fa- Likewise, for models in which the ALFs have a range of 6q values, but all 
else equal, the subset of ALFs with Oq values within about an order of magnitude of the 
highest ^0 value dominate the constraints, while the contribution of the others is negligible. 
Therefore, when comparing actual multiple- ALF models to the constraints and JF values we 
plot below, one should take ualp to be the effective number of ALFs (the number expected 
to have 6q and A^^p values comparable with the largest value) to compare to the results 
below. 

In Figure [21 we show exclusion regions from dark matter density and CMB isocurvature 
constraints and contours of JF for ualf = 1, 10 and 100 ALFs. In each plot, the ALFs have 
A ALF = 78 MeV, and all ALFs are assumed to have the same misalignment angle ^o- Keeping 
all else fixed, the probability of obtaining a certain set of {6q} values is approximately 
given by H^O' t>ut the observational constraints depend on X](^o)^- Therefore, a set of 
misalignment angles {9^} = {1, 10^^°} is as likely to occur as {9q} = {10^^, 10"^}, but the 
latter is much less constrained. Thus models with 6q = 9^ for all z,j as shown in Figure [2] 
represent the minimum tuning for a given constraint level. In the white region in the lower 
left of the plots, the constraints allow each 9q to be of order 1, so the contribution to T 
comes only from the tuning of e (i.e., the inflationary model). This may be read off the 
logio ^ ^^is included along the top of the plots. 

We see that for a single ALF (Figure [21 top panel), the least fine tuning (the largest 
value of JF) that can be achieved without violating cosmological constraints is ~ 10"^^, at 
low values of the symmetry breaking scale Jalf- This corresponds to keeping 9o near its 
"natural" value of order 1 and tuning Hj to be very low (~ 10^ GeV, corresponding to 
e ~ 10"^^). Alternatively, allowing for tuning of the misalignment angle can be traded for 
the tuning of e (e.g., if 9q = 10~^ then Hj can be ~ 10^^ GeV). 

For multiple ALFs, the white region in which the values of 9q can be 0{1) without 
violating constraints becomes slightly smaller, while over the rest of the plot area, the 
parameter J-' is constrained to exponentially small values. This is because, for multiple 
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ALFs, each misalignment angle must be tuned separately to very small values to avoid 
overproducing isocurvature modes or dark matter particles. For the 10-ALF and 100-ALF 
cases (Figure m middle and bottom panels), the level of tuning needed is many tens of orders 
of magnitude worse than the tuning required to solve the strong-CP problem without an 
axion. If Tualf the figure of merit for riALF axions, the degree of tuning scales roughly as 

Bayesian Analysis 

The exponential fine-tuning required to satisfy observational constraints, as parameter- 
ized by JF, suggests that the axion explanation of the strong CP problem, inflation, and 
string theory are not mutually compatible in the absence of some added selection principle. 
To quantify the incompatibility, we introduce a Bayesian model comparison, which incorpo- 
rates both the tuning measure T for ALFs and the prior probabilities of the paradigms. A 
model comparison relies on the calculation of the Bayes factor, which is defined as 

P(D|Mi)P(Mi) 



B 



Ml, Mo 



(9) 



P{D\Mo)P{Mo)' 

where P{D\Mi) is the posterior probability of the data or observation D in the context of 
a model Mj and P{Mi) is the prior probability of the model. A Bayes factor of 3 to 20 is 
considered to be "positive" evidence against model Mq, 20 to 150 "strong" evidence against, 
and greater "very strong" 32|. In the cases of interest here, Bmi,Mo be exponentially 
large. We use the observed upper limit on the CMB isocurvature mode as the data Diso- 
The models we compare are: 

• STwALF: String theory with one or more light ALFs, assuming a high-/a QCD axion 
and inflation, based on the notion that it is difficult to produce a QCD axion without 
also producing light ALFs. 

• no STwALF: No string theory or string theory with no light ALFs other than a 
high- /a QCD axion, and inflation. 

• infi: Inflation, assuming string theory (with ALFs) and a high-/a QCD axion. 



no infi: An alternative to 
mum (e.g., the cyclic model 



inflation in which the axion is never excited from its mini- 
33|), assuming string theory and a high-/a QCD axion. 



cixion: The high-/a QCD axion, assuming inflation and string theory (with ALFs). 
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FIG. 2: Minimum degree of fine tuning required to satisfy cosmological constraints (as measured 
by the figure of merit J-) for (top) 1, (middle) 10, and (bottom) 100 ALFs with instanton scales 
of A ~ 100 MeV. The blue hatched region labelled "overdense" indicates the region excluded by 
dark matter density measurements for any value of the misalignment angle ^o- For each figure we 
also plot the region that is ruled out for all values of by the CMB isocurvature fraction (red). In 
the white region to the lower left, ^ 1 (i-e., no tuning of the misalignment angle is necessary). 
In the middle and bottom panels, for the purpose of illustration, we assume that each ALF has 
roughly the same value of ^o- 

• no cixion: No QCD axion (and an unsolved strong-CP problem), assuming inflation 
and string theory (with ALFs). 

We have shown that for a single axion or ALF, the probability of the model agreeing with 
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the data is T . For the purpose of the Bayesian model comparisons, we are considering the 
combined effect of ALFs and the QCD axion, so we must also factor in the QCD axion's 
misalignment angle, 9^'"^ . Therefore, 

P(A.o I STwALF + axion + infl) = x (10) 

where 6'^'"^ is the value of the QCD misalignment angle required for it, too, to evade the 
isocurvature constraint. Using this, we can calculate Bayes factors expressing the likelihood 
of each of the three paradigms when the other two are assumed. 

For the purpose of illustration, we assume fa ~ 10^^ GeV (close to the string scale) for 
the QCD axion and for ALFs. The value of B decreases somewhat as /„ decreases, but the 
qualitative result is unchanged. A Bayesian model comparison of string theory with light 
ALFs versus an alternative, assuming the existence of the (high- /a) QCD axion and inflation 
yields a Bayes factor of 



P(no STwALF) P{Diso 


no STwALF) 


P(STwALF) P(A.o 


STwALF) 



-B(no STwALF, STwALF) |infl, axions — — n/crr.,, A t t7\ TTTF, lorp... a t t7\ — ("'-■'-) 

P(no STwALF) e^^^O'^^^ 



P(STwALF) 
P(no STwALF) IQ-^^ 



P(STwALF) 10-7(10-4)"ALF 



(12) 
(13) 



^ P(no STwALF) , 
^ P(STwALF) ^ ) 

where e^^^ is the epsilon value required for the QCD axion alone to evade the isocurvature 
constraint and 10~^(10~'^)""*^^ is an approximate expression for the maximum value of JF. As 
for the ratio of model priors, P(no STwALF) /P(STwALF), the current view is that string 
theory is highly favored compared to alternatives and that string theory with ALFs is more 
likely than string theory without. Therefore, this ratio is expected to be less than one, but 
not sufficiently small to overwhelm the exponential factor (10^)"-*^-^. It is one thing to say 
the theory is strongly favored, but it is another to say the chance of an alternative is one in 
ten thousand or exponentially worse. In the equation above, we also have P(i5jso|no ST) = 
(QCDqQCD ^ since, without string theory, ALFs are not produced, so the probability 



is the JF value we would calculate for a single QCD axion 29|. 

By a similar analysis, we find 

P(no infi) P{Diso\no infi) ^ P(no infi) 
P(infi) P(A.o|infi) ~ P(infi) 



-D(no mfi,mfi)|STwALF, axions — — T-,r:^a\ TTTT^ \:^a\ ^ — {^^) ' Ut'j 
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where we take P(iI'jso|no infl) = 1 because neither 60 nor e are constrained by the CMB 
isocurvature mode hmit; and 

P(no axion) 1 -^2/ 4\n 

B(no axion,axion)|infl, STwALF = p^-g^^.^^^^ X ^ ^ ^QCD ~ ) " (l^j 

In this case, we have chosen a model prior that maximally favors the QCD axion as a solu- 
tion to the strong-CP tuning problem, P(no axion)/P(axion) ^ lO"^*^, since no attractive 
alternative exists at present. We have taken P(Djso| axion) ~ J^x^J^^, where 6^"^^ < 10~^ 
is the tuning of the QCD axion misalignment required to evade the isocurvature constraint. 
All three Bayes factors, obtained by making conservative estimates for each contribution, 
become exponentially increasing as ualf increases. 
Discussion 

We have considered the fact that string theory models incorporating a QCD axion typi- 
cally produce additional ALFs with Jalf ~ 10^^ GeV and masses similar to that of the QCD 
axion. Although string theory may allow some vacua with no additional light ALFs, they are 
atypical and not favored by any microphysical arguments. We have shown that the multi- 
ALF models inevitably require a fine tuning of misalignment angles and the inflationary 
slow-roll parameter e that is many orders of magnitude worse than that faced by the QCD 
axion alone, as measured by the figure of merit (JF). In fact, since the ALF misalignment 
angles must each be independently tuned to very low values simultaneously to evade cos- 
mological constraints, the tuning is roughly exponential in the number of additional fields. 
In addition to the tuning in the ALF values of Oq, the effect of inflationary perturbations on 
ALF particle production requires that Hj and hence e be tuned to a much smaller value than 
is necessary to solve the cosmological flatness and horizon problems. Our Bayesian analysis 
quantifies the conclusion that, absent some strong selection principle, inflation, the QCD 
axion explanation for the strong-CP problem, and string theory (as currently understood) 
are not mutually compatible. 

Anthropic reasoning has been invoked in the past as a possible selection principle. How- 



ever, as emphasized in (29|] for the case of the QCD axion alone (no additional ALFs), the 
anthropic principle does not alleviate the problem because most of the anthropically allowed 
range of parameters is ruled out observationally by the isocurvature constraint. That is, 
even after restricting parameters to the range compatible with habitability (based on the 
density constraint), additional tuning is then required to prevent the generation of isocur- 
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vature perturbations whose amplitudes are inconsistent with current measurements of the 
CMB but consistent with habitabihty. This additional tuning is of a magnitude comparable 
to that of the strong-CP problem the axion was invented to solve. 

The same applies to the case of multiple ALFs. For example, for ~ 10^^ GeV, as 
suggested by string theory, the maximum value of T in the parameter space that is anthrop- 
ically allowed but ruled out by the isocurvature constraint is much larger than the maximum 
T in the observationally allowed space. Moreover, in some versions of the anthropic princi- 
ple, it is argued that any anthropically selected variables (such as the dark matter density) 
should have values near the maximal consistent with human existence, also known as the 



anthropic boundary 
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351 ]. If the axionic dark matter density is taken to be near the 



anthropic boundary, large isocurvature amplitudes are strongly preferred, which is at odds 
with observation. 

Since anthropic reasoning does not force the isocurvature mode fraction or the energy 
scale of inflation to be in the narrow range compatible with observations, a compelling 
argument from fundamental physics is needed. Alternatively, at least one of the three ideas 
- the axion solution to the strong CP problem, inflationary cosmology, or string theory - 
must be abandoned. 
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